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ABSTRACT 

Progress  is  reported  in  the  development  and  implementation  of  a  spectral  method 
for  constructing  the  adiabatic  electronic  potential  energy  surfaces  of  large  het¬ 
erogeneous  aggregates  of  physically  or  chemically  interacting  atoms.  The  work 
carries  forward  a  program  of  study  initiated  some  years  ago  by  W.  MofEtt  to 
describe  aggregate  electronic  structure  solely  in  terms  of  the  properties  of  atoms 
and  their  pairwise  interactions.  Necessary  and  sufficient  conditions  are  described 
for  construction  of  the  transformation  from  diatomic  to  spectral-product  states, 
and  of  the  corresponding  pair-interaction  Hamiltonian  matrices  required  in  the 
development,  employing  representatives  of  the  diatomic  spectral  response  op¬ 
erator  obtained  from  Fourier  expansion  of  the  Coulombic  interactions.  These 
conditions  help  to  advance  the  development  of  the  spectral  method  as  a  prac¬ 
tical  computational  tool  for  electronic  structure  determinations  employing  con¬ 
ventional  irreducible-symmetry  methodology.  The  approach  is  illustrated  with 
computational  studies  of  avoided  crossings  in  the  potential  energy  surfaces  of 
small  van  der  Waals  bonded  NzAtn  aggregates,  which  serve  as  useful  prototypes 
for  understanding  selected  attributes  of  cryogenic  high  energy  density  matter 
(HEDM).  An  outline  is  provided  of  work  in  progress  on  general  implementation 
of  the  spectral  method  and  of  its  application  to  cryogenic  fuels  and  oxidizers 
seeded  with  trace  metals  which  may  exhibit  improved  combustion  performance. 


®  Work  supported  in  part  a  grant  from  the  U.S.  Air  Force  Office  of  Scientific  Research  to 
Indiana  University. 

AFOSR  University  Resident  Research  Professor,  1996-98. 
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I.  Introduction 


Development  of  high-speciiic-impulse  cryogenic  fuel-oxidizer  combinations  employing  metal 
or  other  seedant  materials  is  an  important  and  continuing  focus  of  the  USAF  high  energy 
density  matter  (HEDM)  program  [1].  Accurate  adiabatic  electronic  potential  energy  sur¬ 
faces  can  provide  information  useful  for  understanding  the  structure,  stability,  and  other 
attributes  of  such  heterogeneous  systems.  Accordingly,  methods  for  constructing  ab  initio 
potential  energy  surfaces  which  are  applicable  to  large  aggregates  of  physically  or  chemi¬ 
cally  interacting  atoms  are  presently  under  development  employing  a  previously  described 
spectral  method  devised  specifically  for  this  purpose  [2]. 

In  the  present  report,  development  and  implementations  axe  described  of  the  spectral 
method,  and  illustrative  computational  applications  are  presented  in  cases  of  avoided  cross¬ 
ings  in  the  potential  energy  surfaces  of  physically  bonded  (van  der  Waals)  inert-gas  Atn 
aggregates  seeded  with  metal  Na  radicals  [3].  The  general  spectral  method,  which  carries 
forward  a  program  of  study  initiated  some  yeaxs  ago  by  W.  Moffitt  [4]  and  developed  and 
modified  by  others  [5-10],  employs  a  direct  product  of  complete  sets  of  atomic  spectral 
eigenfunctions  to  represent  the  total  adiabatic  electronic  wave  functions.  The  representa¬ 
tion  of  the  electronic  degrees  of  freedom  obtained  in  this  way  is  not  symmetrical  [11],  and 
individual  terms  in  the  many-electron  product  basis  axe  not  exphcitly  antisymmetric  in  the 
electronic  spin  and  spatial  coordinates  [12,13].  The  Hamiltonian  matrix  obtained  in  the 
spectral-product  representation  is  rigorously  pairwise  additive  in  the  individual  atomic- 
interaction  matrices,  which  can  be  obtained  by  performing  calculations  on  a  hmited  num¬ 
ber  of  appropriate  diatomic  molecules  employing  conventional  computational  methodology. 
This  attribute  of  the  spectral  method  can  provide  significant  computational  advantages  for 
laxge  aggregates  of  interacting  atoms,  although  issues  associated  with  convergence  to  an¬ 
tisymmetric  states  must  addressed  in  the  course  of  its  general  implementation. 

Methods  are  described  here  for  constructing  the  required  spectral-product  basis  and  the 
associated  pairwise- additive  atomic-interaction  Hamiltonian  matrices  from  diatomic  molec- 
ulax  eigenstates  and  previously  defined  spectral  response  matrices  [2]  calculated  employing 
conventional  computational  techniques  at  finite  and  asymptotic  (i?  — y  oo)  internucleax 
separations.  These  methods  provide  an  explicit  prescription  for  constructing  the  unitary 
transformation  from  diatomic  to  spectral-product  representations,  advancing  the  develop¬ 
ment  of  the  spectral  method  as  a  computational  tool  for  electronic  structure  determina¬ 
tions.  The  attributes  of  the  approach  are  illustrated  quantitatively  with  calculations  of 
avoided  crossings  in  the  potential  energy  surfaces  of  JVaArjv  clusters  employing  diatomic 
states  constructed  for  this  purpose  in  a  finite  representational  basis. 

Aspects  of  the  formal  theoretical  development  are  briefiy  reviewed  in  Section  II,  the  neces¬ 
sary  and  sufficient  conditions  for  explicit  construction  of  the  unitaxy  transformation  from 
diatomic  to  spectral-product  basis  representations  axe  reported  in  Section  III,  and  illustra¬ 
tive  applications  to  avoided  crossing  situations  in  (JVaArj\f)  van  der  Waals  aggregates  are 
reported  in  Section  IV.  Concluding  remarks  and  an  outhne  of  work  in  progress  on  general 
implementation  of  the  spectral  method  and  its  application  to  physically  and  chemically 
bonded  HEDM  systems  is  given  in  Section  V. 
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II.  Review  of  the  Theoretical  Development 

The  theoretical  development  employs  a  direct  product  basis  of  the  spectral  eigenstates  of 
the  individual  atoms  to  represent  the  correctly  antisymmetric  aggregate  states,  which  basis 
provides  certain  advantages  in  the  study  of  interacting  atoms  [2].  Specifically,  when  the 


product  basis 

$(i;2; . . .  ;n)  =  ®  #^^^(2)  0  •  •  •  (1) 

is  employed  in  a  conventional  linear  variational  solution  of  the  many-electron  Schrodinger 
equation  for  an  aggregate  of  atoms  1  to  iV,  the  Hamiltonian  matrix  takes  the  form 

H(R)  =  {  H<”'  +  '£{13  >  a)  ,  (2) 

Q!=l  ^  /3=1 

where 

H(“)  =  |e(“)0I(^)0I(^)0---I(^)}^,  (3a) 

are  the  atomic  terms,  and 

0  0)  0  •  •  >  (4a) 


are  the  pair-interaction  terms.  In  Eqs.  (1)  to  (4),  is  the  atomic  Hamiltonian 

operator  for  the  atom  a  containing  the  set  of  electrons  i,  y(“’^)(i;j)  is  the  Coulombic 
interaction  operator  for  atoms  a  and  (3,  #(“’^)(*;  j)  is  the  spectral-product  basis  for  the 
indicated  pair,  a  semi-colon  is  employed  to  separate  sets  of  distinguishable  electrons  (i;y), 
the  bracket  symbol  {•  •  •}q  imphes  choice  of  a  particular  ordering  rule  for  the  indices  of  the 
direct  product  (0)  functions,  and  the  vector  R  specifies  the  atomic  positions. 

Because  the  representation  of  Eqs.  (1)  to  (4)  is  not  symmetrical  [11],  the  individual 
products  in  the  basis  of  Eq.  (1)  are  not  totally  antisymmetric  [12,13],  and  the  Hamiltonian 
matrix  of  Eq.(2)  can  contain  terms  which  do  not  correspond  to  physically  meaningful 
(totally  antisymmetric)  aggregate  eigenstates.  These  can  be  removed  by  performing  a 
series  of  pairwise  unitaxy  transformations  on  the  spectral-product  basis  and  the  associated 
Hamiltonian  matrix  which  leave  the  physical  eigenvalues  and  vectors  invariant,  but  which 
allow  unphysical  components  to  be  discarded.  The  result  of  these  transformations  is  to 
replace  the  matrix  elements  of  Eq.  (4b)  with  the  generally  different  quantities  vip’^^(J2a/3) 
given  by  the  expression 

v(“-^)(H,/3)  =  T(“-^)(Jla^)t  •  •  T(“-^)(^a/3),  (5a) 


84 


where 

=  {«<“■«(<, (56) 

is  the  pair- interaction  matrix  in  the  physically  meaningful  portion  of  the  diatomic  eigen- 
state  basis Here,  where 

is  a  product  of  Wigner  rotation  matrices  for  the  atoms  a  and  /?,  the  angles  Woi/3  specify 
the  direction  of  atom  13  from  atom  a,  and  is  the  unitary  transformation  from 

the  physically  meaningful  (antisymmetric)  diatomic  states  to  the  spectral-product  basis  in 
a  pair  of  rotated  atomic  coordinate  systems  having  co-lineax  quantization  axes  [7]. 

III.  Construction  of  the  Pairwise  Rotations 

Necessary  and  sufficient  conditions  for  ab  initio  construction  of  the  pair  transformations 
are  obtained  from  the  expression 


(6a) 

for  the  interaction  matrix  where 


P 

is  the  previously  defined  spectral  response  matrix  in  the  diatomic  basis  [2],  and  the  sum 
over  particles  (p)  in  Eq.  (6b)  includes  all  electrons  {qp  =  — e)  and  nuclei  {qp  =  +Zae  and 
+Zpe)  in  the  molectde.  Equations  (6a)  and  (6b)  follow  from  the  Fourier  transform  of  the 
interaction  operator  appearing  in  Eq.  (5b).  In  the  hmit  of  large  internuclear  separation 
{Rafi  OO) 


ql,“'^>(k.  ^  oo)  ^  M*  .  q<“’«(k,  •  M.  (7) 

where  qlp'^'(k,  is  the  spectral  response  matrix  in  the  pair  product  basis,  and  M  is  a 
known  constant  matrix  that  accounts  for  possible  diatomic  state  degeneracy  in  the  laxge 
Rafi  hmit.  Accordingly,  the  calculated  diatomic  eigenstates  axe  seen  to  provide  the  spectral 
response  matrix  in  both  diatomic  and  spectral-product  representations,  which  matrices  can 
be  employed  to  determine  the  unitary  transformation  matrix  U^“’^^(Raj3)  that  connects 
them 


q<“’«(k,H.p)  =  ■  q^“'«(k,ij.^)  •  (8) 

It  is  seen  from  the  expression  of  Eq.  (6a)  that  the  transformation  of  Eq.  (8)  also  provides 
the  transformation  of  Eqs.  (5).  Because  the  transformation  matrix  U^“’^)(i?a/j)  satisfying 
Eq.  (5a)  also  satisfies  Eq.  (8),  the  latter  provides  necessary  and  sufficient  conditions  for 
determination  of  the  transformation  from  the  physically  meaningful  diatomic  states  to  the 
spectral-product  basis. 
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IV.  Illustrative  Applications  -  NaAr^  Clusters 

Applications  of  the  foregoing  development  are  reported  for  physically  bonded  Ar^  clusters 
seeded  with  JVa  atoms.  For  this  system,  the  rotations  [Eq.  (8)]  axe  solved  in  the  k  — >•  0 
limit,  appropriate  for  long-range  interactions,  in  which  case  Ra^)  becomes  the 

diatomic  dipole  transition  matrix  ^id{Rap)i  a-nd  qlp’^^(k.  Rap)  becomes  the  corresponding 
dipole  transition  matrix,  fisp{Rap),  in  the  spectral-product  basis.  Solution  of  Eq.  (8)  is 
obtained  in  the  form 


where  is  the  unitary  matrix  that  diagonalizes  fJidiRa^)  and  vip‘^\Ra^)  is 

the  uiiitary  matrix  that  diagonahzes  fisp{Roci3)-  Equation  (9)  is  valid  in  the  hmit  that 
closure  has  been  achieved  in  the  diatomic  and  spectral-product  representations  and  the 
eigenvalues  of  /x<i(i2a^)  ^ind  fispiRa/d)  3*re  accordingly  identical. 

Employing  the  transformation  of  Eq.  (9)  in  Eq.  (5a)  gives  the  final  expression 

)'  •  (10) 

for  the  pair-interaction  matrices,  where  the  substitution  — 

has  been  employed,  'E^^'^\Rap)  is  the  diagonal  matrix  of  diatomic 
energies,  and  is  its  Rap  oo  limit. 

In  Fig.  1  axe  shown  diatomic  potential  energy  curves  for  the  ground  and  low-lying  excited 
states  of  NaAr  calculated  as  functions  of  internuclear  distance  employing  multi-reference 
configuration-interaction  methods  [14].  Evidently,  the  excited  F  and  H  states  undergo  an 
avoided  crossing,  whereas  the  other  states  within  a  symmetry  class  (S,  11,  A)  are  generally 
well  sepaxated  at  aU  Rap-  The  associated  electric  dipole  and  transition  matrix  elements 
required  in  Eq.  (8),  shown  in  Figs.  2  and  3  for  selected  S  and  11  states,  are  seen  to  be 
strong  functions  of  internuclear  separation,  consequent  of  atomic-product-state  mixings  at 
aU  Rap-  Also  shown  in  these  figures  axe  the  eigenvalues  of  the  matrices  Hd{Rap)  for  S 
and  n  states,  which  are  significantly  less  dependent  upon  internuclear  separation  than  are 
the  matrix  elements  themselves.  Moreover,  these  eigenvalues  axe  in  very  good  accord  with 
those  of  the  matrices  fj,sp{Rap),  which  are  also  shown  in  the  figures,  and  which  should  be 
entirely  independent  of  Rap-  In  the  limit  of  closure  in  the  diatomic  and  spectral-product 
basis  sets,  the  eigenvalues  of  the  two  matrix  representatives  of  the  transition  moment 
operator  become  identical  and  independent  of  Rap- 

The  transformation  matrices  obtained  for  Na-Ar  interactions  from  the  results  of  Figs.  1 
to  3  and  Eqs.  (6)  to  (10)  with  k  0  are  employed  in  calculations  of  potential  energy 
surfaces  for  the  aggregates  of  iaterest,  whereas  the  Ar-Ar  interactions  are  treated  as  pair¬ 
wise  additive  in  the  ground-state  potentials.  Note  that  the  matrix  M  becomes  the  unit 
matrix  for  the  states  of  interest  in  NaAr.  In  Fig.  4  are  shown  sections  of  the  potential 
energy  surfaces  of  iVaAr2  {G2v  symmetry)  as  functions  of  the  Ar-Na-Ar  angle  for  states 
associated  with  Na  3s  -t  3d,  4p  excitations.  Also  shown  in  the  figure  are  the  corresponding 
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states  obtained  from  the  present  development  employing  a  unit  matrix  for  the  required 
transformation  [U(“’^)(i2a/3)  =  1],  which  corresponds  to  the  so-called  Balling  &  Wright 
approximation  commohly  employed  in  coUisional  broadening  studies  and  trapped  radical 
spectroscopy  [3].  Evidently,  the  spectral  method  provides  avoided  crossings  in  the  calcu¬ 
lated  ^Ai  and  ^£2  energy  surfaces  associated  with  JVa  3a  3d,  4p  excitations,  whereas  the 
simpler  model  does  not.  Similar  observations  apply  to  the  more  comphcated  energy  sur¬ 
faces  for  larger  JVaArjv  aggregates,  also  obtained  from  the  present  development  employing 
the  single  Na-Ar  transformation  matrix  of  Eqs.  (9)  and  (10).  These  results  suggest  that 
non-additive  contributions  to  the  potential  energy  surfaces  of  large  van  der  Waals  aggre¬ 
gates  can  be  determined  quite  generally  in  terms  of  the  underlying  pairwise-interaction 
matrices  by  employing  the  computational  prescription  of  Eqs.  (9)  to  (10)  for  constructing 
the  required  pair  transformation  matrices.  In  the  case  of  chemical  interactions,  the  more 
general  spectral  response  matrices  of  Eqs.  (6)  to  (8)  are  preferred  over  the  dipole  matrices 
employed  for  van  der  Waals  aggregates. 

V.  Concluding  Remarks  -  General  Implementation  and  HEDM  Applications 

The  spectral  method  and  procedures  described  here  for  computational  implementations 
have  been  employed  in  studies  of  the  effects  of  non- additivity  in  the  potential  energy 
surfaces  on  the  structure  and  optical  absorption  spectra  of  NaAr^  [W=  2  to  147]  van  der 
Waals  aggregates,  results  that  will  be  reported  subsequently.  Some  indication  of  the  nature 
of  the  approach  is  provided  by  the  illustrative  results  reported  in  Figs.  1  to  4,  which  employ 
a  limit  (k  -t  0)  appropriate  for  long-range  interactions.  Evidently,  avoided  crossings  are 
obtained  from  the  spectral  method  in  the  potential  surfaces  of  these  aggregates  employing 
accurate  diatomic  potential  curves. 

When  the  development  of  Eqs.  (6)  to  (10)  is  employed,  rather  than  the  dipole  hmit  em¬ 
ployed  in  the  apphcations  reported  here,  the  spectral  method  and  its  implementation  is 
generally  applicable  to  both  physical  and  chemical  interactions.  Evaluation  of  the  spectral 
response  matrix  of  Eqs.  (6a)  and  (6b)  in  the  large  electronic  basis  sets  required  for  di¬ 
atomic  calculations  will  entail  development  of  appropriate  computer  code  routines  for  this 
purpose.  Additionally,  adoption  of  previously  described  Stieltjes  methods  for  construction 
of  complete  sets  of  spectral  states  will  allow  apphcations  of  the  spectral  method  in  the  hmit 
of  closure  [15].  These  issues  are  presently  under  study  and  wih  be  topics  of  subsequent 
HEDM  reports. 

Apphcations  of  the  spectral  method  axe  in  progress  to  potential  cryogenic  HEDM  systems 
which  may  improve  the  specific  impulse  of  H2  /  O2  propulsion  systems.  Specifically,  the 
potential  energy  surfaces  of  Li(H2)N  and  A1(H2)n  aggregates  are  under  study  in  order  to 
clarify  the  natures  of  the  physical  (van  der  Waals)  trapping  sites  in  these  cases,  and  to 
determine  stabihties  of  the  trapping  configurations  to  chemical  reactions  that  can  produce 
the  stable  products  LiH  and  AJH  when  the  aggregates  are  subjected  to  thermal  or  other 
influences.  Additionally,  in  view  of  the  importance  of  cryogenic  He^  clusters  as  potential 
vehicles  for  pickup  and  transport  of  trace  metals  into  cryogenic  fuel  or  oxidizer  matrices, 
studies  of  the  potential  energy  surfaces  of  LiHejv  and  AJifejv  aggregates  will  be  undertaken. 
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